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摘 要 : 为 实现 在 低 信 骂 比 的 情况 下 对 载 频 进行 参数 估计 ， 提 出 了 一 种 算法 。 首 先 ， 根 据 循环 谱 理 论 计 算出 信号 的 循 
环 谱 ; 然后 ， 采 用 压缩 感知 的 方法 重新 构造 信号 的 循环 谱 切面 ; 最 后 ， 根 据 信 号 循环 谱 切面 中 的 离散 谱 线 位 置 与 信号 
的 载 频 的 关系 采用 平均 法 来 计算 载波 频率 。 数 值 仿真 显示 在 信号 稀疏 度 较 高 的 情况 下 ， 对 载波 频率 的 估计 精度 能 控制 
在 103 Hz 以 下 。 在 低 信 噪 比 的 情况 下 ， 该 算法 能 够 有 效 地 实现 对 信号 载 频 的 估计 ， 在 工程 上 具有 一 定 的 实用 价值 。 
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Signal carrier frequency estimation based on compressed sensing and cyclic Spectrum theory 
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Abstract: Parameters estimation of the modulated signal, such as carrier frequency, Symbol rate, etc. , is the precondition for 
signal modulation recognition. This paper propose the algorithm to realize the parameter estimation of the carrier frequency 
under the condition of low SNR. Firstly, the algorithm calculates the cyclic spectrum of the signal based on the cyclic spectrum 
theory. Then, it rebuilds the cyclic spectrum of the signal by the compressed sensing method. Finally, it calculates the carrier 


frequency using the average method according to the relationship between the position of the discrete spectral line in the signal 


cycle spectrum and the carrier frequency of the signal. The simulation results show that the estimation accuracy of carrier 
frequency can be controlled below 103 Hz under the condition of high signal sparsity. In the low SNR environment, the algorithm 


can effectively estimate the carrier frequency of the signal, which has certain practical value in engineering. 
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使 大 量 信号 呈现 出 循环 平稳 性 ， 因 此 谱 相 关 理 论 在 信号 检测 、 
参数 估计 、 系 统 辨 识 、 调 制 识 别 等 现代 信号 处 理 诸 多 领域 中 得 

近 几 十 年 ， 随 着 现代 通信 技术 的 高 速 发 展 ， 调 制 方式 识别 到 了 越 来 越 深 入 的 研究 。 其 中 ， 载 波 频率 和 符号 速率 在 解 调 中 
和 信号 参数 的 估计 在 民用 或 者 军事 领域 都 扮演 着 非常 重要 的 角 ”起 着 重要 作用 。 在 军事 侦察 、 住 宅 无 线 电 管制 和 软件 接收 机 等 
色 ， 例 如 防御 灾害 、 地 质 、 海 洋 的 探查 、 电 子 对 抗 ， 情 报 侦察 ”” 非 合作 通信 环境 中 ,这 两 个 参数 是 解码 信号 的 必要 条 件 . 因 此 ， 
等 领域 等 。 而 信号 调制 参数 的 估计 ， 如 调制 识别 ， 是 成 功 识 别 ”在 过 去 几 十 年 的 研究 中 ， 越 来 越 受到 人 们 的 重视 。 


0 引言 


号 调制 方式 的 重要 前 提 ， 通 常情 况 下 ， 信 号 调制 参数 包括 : 估计 载波 频率 的 方法 有 很 多 ， 在 基于 频 域 分 析 方 面 ， 文 献 
载波 频率 、 码 元 速 言 号 带宽 等 。 [2] 的 周期 图 法 和 文献 [3] 的 频率 居中 算法 适用 于 具有 很 强 的 载 
谱 相关 理论 由 Gardner 等 人 一 深入 研究 并 发 展 ， 其 优点 主 。 波 功率 和 对 称 的 功率 谱 密 度 信号 。 在 基于 时 域 分 析 方 面 ， 文 献 


要 表现 在 分 辩 率 高 、 抗 干扰 能 力 强 、 双 频率 域 信 息 丰 富 且 易于 [3] 过 零点 算法 ， 文 献 [4] 带 状 谱 估计 算法 (SSCA) 计算 简单 、 
提取 等 ， 适 用 于 循环 平稳 信号 的 分 析 。 在 通信 、 雷 达 、 遥 测 、 计算 量 小 ， 但 是 对 噪声 敏感 。 在 基于 小 波 方面 ， 文 献 [5] 和 文献 
水 声 、 电 子 对 抗 等 领域 中 ， 由 于 某 些 人 为 的 或 无 意 的 调制 ， 而 [6] 介 绍 了 关于 小 波 变换 以 及 基于 负 科 最 大 化 的 算法 ， 去 噪 性 能 
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好 ， 但 计算 复杂 ， 计 算 量 大 。 文 献 [7] 提 出 了 计算 循环 谱 和 截面 进 


行 降 维 处 理 ， 使 得 计算 量 大 大 降低 ， 但 是 适应 性 差 。 在 改善 特 
环 谱 算法 方面 , 文献 [8, 9] 通 过 分 析 循 环 谱 矩 阵 ， 一 方面 改善 特 
环 谱 和 矩阵 的 计算 方法 ， 更 方便 的 提取 特征 值 ， 男 一 方面 对 循环 


谱 进行 降 维 处 理 ， 使 得 循环 谱 计 算 量 大 大 降低 。 文 献 [10~12] 提 
出 对 信号 采样 值 进行 数据 分 段 ， 而 且 每 段 数据 与 前 段 数 据 重生 
一 部 分 ， 得 到 每 段 信号 的 循环 谱 ， 再 进行 平均 处 理 ， 得 到 信号 
的 循环 谱 ， 使 得 噪声 影响 降低 ， 提 高 了 信 噪 比 门限 。 对 于 在 有 
限 数据 情况 下 ， 改 善 了 算法 的 有 效 性 。 但 是 其 分 段 次 数 需 
操纵 , 使 得 算法 稳定 性 不 高 , 且 增 大 了 算法 的 计算 量 。 文 献 [13， 
14] 提 出 了 基于 高 阶 累积 量 与 循环 谱 相 结合 对 载 频 进行 精确 佑 
计 的 方法 ,该 方法 利用 了 高 阶 累积 量 在 抗 噪声 干扰 方面 的 优势 ， 
提取 信号 的 特征 量 。 算 法 精确 度 高 ， 抗 噪 强 。 但 计算 量 较 大 。 


在 信号 处 理 起 至 
言 号 调制 分 类 识别 的 前 
率 越 来 越 高 ， 


谢 瑶 ,， 等: 


5 


I 中 流 碟 柱 的 作用 。 
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(D=Re(zr)e ”ar 


其 中 ，S8( 门 被 称 为 循环 谱 密度 函数 。 


压缩 感知 
奈奈 斯 特 采 样 定理 作为 模拟 信号 和 数字 信号 的 桥梁 ， 一 直 


存储 、 传 输 产生 巨大 上 


的 有 


提出 压缩 感知 理论 (c 


盲 号 的 采用 处 天 
和 i 提 ， 随 着 奈奈 斯 特 信号 采样 要 求 的 采样 


(5) 


是 实现 数字 


且 信 号 离散 后 会 生成 大 量 的 离散 数据 ， 给 数据 
E 力 。Candes，Donoho 和 Tao 等 人 呈 间 
ompressed sensing,CS )， 为 这 一 问题 的 解 


多 化 算法 将 信号 进 


于 循环 谱 在 频 域 上 表现 出 稀疏 性 ， 文 献 [15] 提 出 了 采用 压缩 
感知 的 方法 对 循环 谱 进 行 重 构 ， 消 除 循 环 谱 中 的 元 余 信息 ， 提 
取信 号 中 的 有 效 信息 。 但 是 需要 较 多 的 测量 次 数 ， 且 重建 精度 
相对 较 低 。 
针对 现 有 算法 存在 的 这 些 问题 ， 本 文 提出 了 一 种 基于 分 段 
正 交 匹配 追踪 (Stagewise OMP) 的 算法 。 在 算法 的 欠 代 过 程 中 ， 
每 次 可 以 选择 多 个 元 素 ， 大 大 减少 了 信号 重 构 的 时 间 。 此 外 ， 


pp 


算法 采用 门限 值 的 方式 来 确定 元 素 ， 少 了 元 素 个 数 的 限定 ， 避 


行 精确 


决 提供 了 新 的 思路 ， 即 在 信号 具有 稀 玻 性 的 情况 下 ， 可 以 通过 
打破 了 奈 奎 斯 特定 理 的 限制 ， 


Im 
Ia 
起 


有 效 的 降低 了 采样 速率 。 


于 超 宽 乔 


采样 设备 的 依赖 ""。 当 信号 在 某 一 域 是 
稀 琉 的 ， 贝 


和 信号 而 言 ， 可 以 通过 寻 


找 

的 压缩 数据 实现 对 信号 无 失真 的 重 构 恢 复 ， 解 除了 信和 号 对 前 端 
稀 
的 


上 可 以 通过 一 个 满足 一 定 条 件 


ij 且 投 影 后 的 


免 导致 元 素 信息 的 见 余 和 缺失 ; 本 算法 不 依赖 于 信号 的 稀 疏 度 ， 
对 信号 的 重 构 有 着 独到 的 优势 。 
1 ”循环 谱 的 计算 方法 

循环 谱 是 信号 平稳 分 析 的 常用 方法 ， 利 用 不 同 的 信号 循 : 
谱 的 差异 可 以 对 参数 不 同 的 信号 进行 识别 。 由 于 平稳 噪声 不 具 
有 循环 平稳 特性 ， 因 此 循环 谱 分 析 可 以 在 很 大 程度 上 将 信和 号 与 
噪声 区 分 开 ， 在 信号 分 析 中 具有 很 大 的 优越 性 “ ”。 

设 x() 的 一 阶 统 计 特 性 均值 和 二 阶 统计 特性 自 相关 函数 具 
有 某 些 周期 工 的 特性 ， 则 称 之 为 广义 周期 平稳 过 程 。 均 值 及 自 
相关 函数 分 别 为 ”“: 


A 


mAD=m.(t+7) (D) 


RtT2,t-7D= RTHT2,t+T -7/2) @Q) 
其 中 : T= K/a，K 为 整数 ，a 为 周期 频率 。 将 周期 函数 
Rx(t,T) 展 开 傅 里 叶 序 列表 达 式 为 


i2xat 


R=2 RWe G) 


aeD 
其 中 : R&(7) 为 循环 自 相关 函数 , 是 循环 频率 ag 和 时 间 延 迟 t 的 函 
数 。 其 表达 式 为 


-i2xat 


RO=limz) Rr/ Ye di (4) 


若 g=0， 则 R%(7) 为 平稳 信号 的 自 相关 函数 。 循 环 相关 函数 
Rg(T) 的 傅 里 叶 变 换 表 达 式 为 


采样 频率 与 信号 


的 采样 频率 是 不 相关 的 ， 也 生 


其 稀 玻 域 ， 利 用 少量 


琉 的 或 者 其 本 身 就 是 
4 测量 矩阵 使 得 投影 后 


号 具有 稀 玻 性 ，T 言 写 包 含 了 原 信 和 号 足够 多 的 
J 以 通过 算法 重 构 恢 复出 原 信 号 。 因 此 ， 压 缩 感知 理论 


1 原始 信号 的 


频率 不 相关 ， 只 是 与 信号 本 身 的 性 质 相 关 。 出 


上 外， 信号 的 


处 理 可 以 不 通过 奈 奎 斯 特定 理 来 采样 信号 ， 


投影 到 一 个 低 维 空间 中 "。 
言 号 的 稀 玻 表示 
的 稀 疏 性 也 可 称 为 信号 的 可 压缩 性 ， 
下 ， 其 表示 的 系数 只 有 少数 不 为 零 ， 
。 信 和 号 的 稀 玻 分 析 方 法 是 在 一 定 的 先 验 条 


据 先 构造 出 的 基 函 数 来 


压缩 感知 


则 x 是 可 压缩 的 。 信 和 号 的 稀 杏 


首先 就 是 要 找到 信和 号 的 稀 玻 域 。 对 于 信号 xE R™， 
E 交 基 的 线性 组 合 来 表示 ， 旧 


x 


1 


rz=> 00D =Y@ (6) 


而 是 通过 算 


是 指 信号 在 某 一 


或 者 是 多 数 系数 


牛 的 基础 上 ， 


可 效 地 表示 数据 。 对 于 信和 号 来 说 ， 


可 以 由 一 组 


其 中 : 下 = [gi1,92,…,9w] 为 正 交 基 和 矩阵 ， 列 向 量 gi 为 基 函 数 ， 


8 = [61,6…,9w]7， 是 系数 向 量 。 如 果 9 只 有 K 个 系数 较 大 ， 


的 信号 在 观测 过 程 中 都 是 非 


如 果 存 在 正 交 矩阵 下 使 得 x 


(7) 


i 性 决定 了 信号 压缩 的 最 佳 效 果 ， 
对 于 信号 x， 在 压缩 处 理 中 只 需要 保留 前 K 个 最 大 的 系数 ， 也 
不 会 导致 信号 的 质量 出 现 严重 下 降 。 
2.2 观测 矩阵 

对 于 信号 xE RN， 
观测 矩阵 e RMXN 对 信号 x 进 
其 表达 式 为 

y= Dx 


是 天 - 黎 玻 ， 


行 观测 ,得 到 ye RN,K<M<N。 


适应 性 的 , 测量 


E 阵 性 质 的 高 低 ， 与 成 功 压缩 信号 和 精确 恢复 重 构 信 号 息 息 相 
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关 。 设计 测量 矩阵 需要 满足 一 定 的 条 件 才能 实现 信号 的 重 构 ， 
2006 年 Candes、Romberg 和 Tao 等 人 提出 了 著 名 的 约束 等 距 
性 (Restricted Isometry Property，RIP) 理 论 ，RIP 准则 限定 了 观 
测 矩 阵 的 设计 ， 成 为 压缩 感知 葛 基 性 理论 ， 为 压缩 感知 理论 的 
准确 性 提供 了 保障 。 即 假定 信号 x 的 投影 稀 疏 向 量 @ 的 长 度 为 
N， 稀 玻 度 为 玉 ， 对 感知 矩阵 A = @ 和 e RYXN， 若 存在 常数 6 E 
(0,1)， 满 足以 下 不 等 式 : 


a-oel. <|A8|. <arale 


则 称 感知 矩阵 A 满足 K 阶 约束 等 距 性 准则 ” ”。 
2.3 ”信号 重 构 
压缩 采样 的 成 功 与 否 ， 关 键 在 于 能 否 成 功 重 构 原 始 信号 。 
在 压缩 感知 理论 中 ， 由 于 观测 算法 的 数量 小 于 未 知 数 的 数量 ， 

信号 X 的 重 构 问 题 是 一 个 欠 定 问题 。 即 ， 输 入 为 观测 值 y 和 感 


三 


2 2 
< 
2 2 


(8) 
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其 中 ,t= 12,…,N, 9$(Tt) = 3 aill(Tt -iT,), ie 


人 多 直 _，，A 为 信号 幅度 ，N 为 样本 数 ，T 为 时 间 间 隔 ，N; 为 


码 元 个 数 ， 罗 为 码 元 宽度 ，[[( ) 为 脉冲 形成 函数 ， 无 为 信号 载 
波 频率 ，68o 为 信号 的 初始 相位 。 假 设 
SE(f) = QreU+a/2)Qre( — a/2) (13) 


其 中 , Qre(f) = 2 人 。 根 据 表达 式 (3)~ (5) 可 以 算出 BPSK 


害 号 的 循环 谱 为 中 
A A A 


Se 2fe(f)ej200] (14) 


其 中 :w = +2f. + kT,，a = kT,，k 是 整数 。 对 于 平稳 白 噪声 ， 
其 谱 密度 集中 在 a=0 区 域 。 在 oz 时 ， 噪 声 的 影响 很 小 。 结 果 


ed 


知 矩 阵 A， 在 y= A8 的 约束 下 求 稀 疏解 或 其 允 近 值 6， 其 方程 
解 有 无 数 个 。 表 达 式 为 : 

6 = arg minllello s.t. y= AQ (9) 
0- 范 数 表示 9 中 非 零 元 素 的 个 数 ， 是 一 个 非 凸 优化 的 问题 ， 


表明 ， 基 于 循环 谱 的 参数 估计 方法 可 以 有 效 地 抑制 噪声 。 

然后 采用 贪 梦 重 构 算法 重新 计算 循环 谱 截 面 52(f)。 贪 梦 算 
法 中 最 常用 的 算法 正 交 匹配 追踪 算法 (OMP), 其 目的 是 通过 在 
多 次 迭代 中 查询 出 最 少 列 向 量 作为 基 向 量 对 原始 输入 的 稀 疏 多 
频带 信号 的 频 域 信号 进行 表示 ， 最 终 通过 这 些 基 向 量 的 加 权 函 


Im| 


方程 难以 直接 求解 。 为 解决 这 一 问题 ， 通 常 利 用 松弛 技术 采用 
凸 优化 问题 对 其 进行 逼近 。 在 压缩 感知 中 ， 会 把 非 凸 的 目标 函 
数 转化 为 凸 函数 。 转 化 后 的 方程 为 

6 = arg minllell s.t.y= A9 


(10) 

言 写 的 重 构 算法 主要 分 为 优化 类 算法 和 贪 禁 算法 。 优 化 类 
算法 包括 基 追 踪 法 (Basic Pursuit,BP )、 山 优化 算法 等 算法 。 贪 
焚 算 法 包括 正 交 匹配 追踪 算法 (COMP)、 正 则 化 正 交 匹配 追踪 算 
法 (ROMP)、 分 段 正 交 匹配 追踪 算法 (STOMP) 等 算法 。 
的 重 构 过 程 就 是 利用 最 优 的 查询 方法 找到 带 有 有 效 频 谱 信 息 
频段 及 其 所 对 应 字典 中 的 位 置 ， 从 而 实现 信号 的 恢复 。 由 于 查 
询 方法 均 具有 较 高 的 复杂 度 ， 具 有 较 高 精度 的 查询 能 力 的 方法 
计算 非常 复杂 ， 需 要 考虑 或 者 附加 的 条 件 很 多 ， 这 就 导致 这 些 


A 


5 村 
东 < 


于 


数 对 信号 进行 重 构 。 算 法 的 过 程 是 在 每 次 迭代 中 查询 所 需要 的 
信息 和 信息 对 应 的 位 置 ， 然 后 重新 残 差 ， 一 直 满 足 迭 代 终 止 条 
件 时 才 会 停止 迭代 。 在 巴 优化 算法 中 虽然 计算 量 大 和 计算 复杂 
度 高 ， 但 是 其 鲁 棒 性 比较 高 ， 而 OMP 算法 虽然 在 计算 速度 上 
相 比 凸 优化 算法 有 着 很 大 的 改善 ， 但 是 其 鲁 棒 性 较 低 。 所 以 在 
结合 凸 优化 算法 和 OMP 算法 的 优点 后 ， 提 出 了 正则 化 正 交 匹 
配 追踪 算法 (ROMP)。 相 比 OMP 算法 ，ROMP 算法 是 在 每 次 
迭代 中 会 任意 选择 出 列 向 量 ， 而 非 只 选择 一 个 列 向 量 。 因 为 在 
随意 选择 的 列 向 量 中 可 能 会 出 现 误 选 的 情况 ， 为 了 避免 每 次 迭 
代 中 过 多 的 误 选 ，ROMP 算法 通过 引入 正则 化 的 方法 对 待 选 的 
矩阵 作 约 束 ， 进 而 降低 误 选 率 、 减 少 迭 代 和 提高 恢复 精度 。 

与 OMP 算法 一 样 ，ROMP 算法 在 信号 重 构 中 也 需要 知道 


算法 在 理论 上 很 占 优势 ， 应 用 到 实际 的 工程 的 优势 不 大 。 在 实 
际 工程 应 用 中 更 倾向 于 贪 禁 算 法 。 贪 禁 算法 计算 复杂 程度 不 高 ， 


但 对 硬件 的 要 求 较 高 。 通 过 尽 可 能 少 的 迭代 次 数 ， 查 询 出 最 具 
代表 性 的 向 量 对 原始 输入 稀疏 多 频带 做 加 权 表 示 ， 然 后 利用 加 
权 表 达 式 对 信号 进行 重 构 ” 。 

3 ” 载 频 估计 

整形 滤波 器 的 特性 ，Q7c( 有 ) 在 全 0 处 取 最 大 值 。 当 a = 0 
时 ， 载 波 频率 可 以 在 其 频 域 被 计算 。 首 先 计算 信号 的 循环 谱 


Q 


0.(f)， 并 求 出 w 二 0 截面 。 


以 MPSK 信号 模型 为 例 ， 假 设 接收 到 的 信号 为 : 
X(t) = S(t) + n(t) (11) 
S(t) 为 MPSK 信号 ，n(b 为 均值 为 0 方差 为 c2 的 高 斯 白 噪 
声 。S(t) 的 数学 表达 式 为 : 
S(t) = Aexp{j2nzf.Tt + 00 + $b(TE)} (12) 


言 号 的 稀 艳 度 来 决定 在 迭代 中 需要 的 最 大 迭代 次 数 。 此 外 ， 
ROMP 算法 也 存在 一 定 的 缺陷 : 当 感 知 输出 矩阵 中 含有 噪声 的 
时 候 ， 会 导致 迭代 无 法 停止 。 在 遇见 这 种 情况 时 ， 就 需要 先 验 
的 稀 朴 度 来 控制 迭代 的 次 数 。 
为 了 进一步 提高 迭代 速率 ， 分 段 正 交 匹 配 追 踪 算 法 
(STOMP ) 是 在 正 交 匹配 追踪 算法 的 思想 基础 上 加 入 了 分 段 查 
询 的 思想 ， 即 在 每 次 迭代 时 设置 一 个 阀 值 ， 当 字典 中 的 列 向 量 
与 残 差 的 相关 度 大 于 该 阔 值 的 列 向 量 均 视 为 正确 的 频谱 支撑 区 
查询 结果 ， 根 据 算 法 残 差 的 相关 度 的 情况 ， 阀 值 的 选择 一 般 在 
2~3"”。STOMP 算法 实现 步骤 如 下 : 
a) 初始 化 : 稀 跑 解 9=0， 感 知 和 矩阵 A = 下， 输出 矩阵 y， 
残 差 向 量 e=y。 支 撑 集 设置 为 空 ， 即 1 = $p， 设 置 闵 值 TH。 
b) 残 差 的 后 向 投影 : wx = 47ek_ 1 根据 字典 的 随机 性 假设 ， 
大 多 数 计算 得 到 的 值 都 是 均值 为 0 的 独立 分 布 的 高 斯 向 量 。 
过 比较 ux 与 阐 值 TH 的 大 小 ， 得 到 相对 应 的 A 的 列 序号 集 
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J。 下 
c) 支撑 集 : 在 迭代 过 程 中 , 每 次 迭代 都 更 新 感知 矩阵 A， m 
4 = Ai_1U qs， 和 支撑 集 I,， 1 = JiUJ。 
d) 更 新 残 差 : 用 最 小 二 乘法 计算 Be = argnynlly - go| 
Ai9il| = (4TA1)-14Ty。 得 到 更 新 后 的 Qi 后， 利用 e: =y 一 y 一 | 
4,9,, 得 到 新 的 残 差 。 
STOMP 和 迭代 次 数 要 小 于 OMP 的 次 数 ， 相 当 于 简化 了 算 2of 
法 ， 而 且 在 每 次 迭代 过 程 中 会 选择 多 个 列 向 量 ， 由 阔 值 的 大 小 a = = 
来 决定 此 次 选择 列 向 量 的 个 数 ,每 次 狗 代 会 选择 出 多 个 基 向 量 ， ew 
收敛 速度 更 快 。 另外， 因为 采用 分 段 的 阐 值 来 选择 基 向 量 ， 所 图 2 全 2 时 重 构 概率 图 
以 信号 的 重 构 误差 也 较 小 。 二 
最 后 ， 搜 索 通 过 贪 禁 算 法 重 构 出 的 循环 谱 截 面 9 里 面 ool 
的 峰值 ， 利 用 频率 平均 法 估算 载 频 。 其 频率 平均 法 公式 如 下 : "| 
二 0 | 


综 上 所 述 ， 载 频 估计 的 完整 算法 流程 如 图 1 所 示 。 
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3 ”人 2.6 时 重 构 概 率 图 
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图 4 t=3 时 重 构 概 率 图 

从 图 2~4 可 以 看 出 ， 当 信号 在 某 一 域 的 稀 疏 性 越 好 ， 需 要 
观测 提取 的 个 数 越 多 ， 信 和 号 被 重 构 出 来 的 概率 就 越 大 ， 当 算法 
的 阀 值 取得 越 大 ， 从 原 信号 中 得 到 的 信息 就 会 越 少 ， 若 需要 保 
证 信号 被 重 构 出 来 的 效果 越 好 ， 则 需要 观测 的 数据 也 就 越 多 ， 
当 阀 值 取 得 小 时 ， 从 原 信号 会 提取 到 宛 余 的 信息 ， 为 保证 重 
构 效 果 ， 需 要 的 观测 数量 也 就 越 多 ; 当 阀 值 为 2.6 时 ， 重 构 效 
果 了 上 略 优 于 立 值 为 2 和 3 的 时 候 。 当 观测 的 数量 达到 一 定 的 数量 
图 1 算法 流程 图 后 ， 信 号 就 可 以 被 完整 的 重 构 出 来 。 

4.2 ”循环 谱 与 重 构 结果 

由 于 平稳 噪声 和 干扰 在 a 0 处 不 呈现 谱 相 关 , 因此 在 利用 
4.1 系数 对 算法 的 影响 a 关 0 处 的 谱 相 关 函 数 对 调制 信号 进行 检测 与 估计 时 , 可 以 完全 
言 号 的 恢复 是 压缩 感知 的 关键 。 信 号 在 某 一 域 表现 的 稀疏 ”摆脱 背景 噪声 的 影响 ， 采 用 BPSK 方式 调制 的 有 噪声 背景 下 ， 
性 越 稀 疏 , 信号 被 重 构 的 概率 越 大 。 对 于 本 文 提 到 的 STOMP 算 在 全 0 的 循环 谱 截 面 上 可 以 很 容易 地 检测 到 2 倍 载 频 , 在 背景 
法 ， 门 限 参数 的 选择 和 观察 数量 的 个 数 对 于 算法 重 构 的 信号 效 ”噪声 下 检测 到 载 频 的 存在 就 完成 了 信和 号 的 检测 。 本 实验 采用 
果 起 着 决定 性 的 作用 ,在 信号 个 数 N=256, 阀 值 取 值 与 稀疏 度 ， BPSK 信号 为 例 ， 其 仿真 参数 为 : 采样 频率 f=1600 kHz, 载 波 频 
观测 个 数 的 关系 如 图 2~ 4 所 示 。 率 太 = 200 kHz， 码 速率 亡 = 10 kbps， 信 噪 比 SNR=-3 dB， 数 


搜寻 重 构图 的 峰值 对 应 的 频率 值 


求 和 平均 得 到 频率 值 
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据 长 度 N=1 000。 使 用 MATLAB 对 BPSK 信号 进行 了 循环 谱 的 
仿真 , 并 对 循环 谱 个 0 和 a = 0 的 截面 仿真 。 其 仿真 图 形 如 图 5~ 
7 所 示 。 


BPSK 循 环 谱 三 维 图 
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图 6 BPSK f= 0 循环 谱 截 面 
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图 7 BPSK a = 0 循环 谱 截 面 

对 于 BPSK 信号 ， 只 有 在 a = 土 2 十 kTb，a = kT。， 的 截 
面 才能 出 现 循环 谱 信号 。 其 中 天 为 载波 频率 ， 思 为 码 元 周期 。 对 
于 实 信 号 ， 在 谱 相 关 图 中 ， 由 于 循环 谱 算法 具有 对 称 性 ， 只 需 
要 知道 a € (0，,+co) 和 fe (0，,+co) 的 部 分 , 就 可 以 知道 循环 
谱 的 图 形 。 

对 循环 谱 和 截面 信号 通过 STOMP 算法 进行 信号 重 构 得 到 的 
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图 形 如 图 8 所 示 
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图 8 BPSK f= 0 重 构 图 
从 图 8 可 以 看 出 ， 由 式 (15) 可 以 得 出 ， 在 f= 0 的 循环 谱 


截面 BPSK 信号 只 在 2 倍 频 才 有 频谱 信号 。 通 过 搜寻 重 构图 内 


最 大 的 两 处 峰值 ， 求 出 相对 应 的 频率 值 ， 就 可 以 得 到 相应 的 载 
频 。 由 上 图 可 以 得 出 测量 的 峰值 为 406KHz 与 -394KHz。 与 循环 
谱 理论 值 存在 6KHz 的 误差 。 由 式 〈16) 得 测 得 到 的 载波 频率 
/各 -200KHz， 算 得 载 频 与 实验 结果 一 致 。 由 实验 结果 
分 析 ， 循 环 谱 对 载 频 的 估计 精度 在 103 Hz。 

5 结束语 


函数 来 研究 非 平稳 信号 ， 被 广泛 地 应 用 于 


检测 算法 中 。 本 


性 ， 


文 分 析 了 应 用 循环 谱 理 论 进行 信号 检测 和 


循环 平稳 理论 主要 是 利用 信号 的 统计 参量 ， 如 均值 和 相关 


参数 提取 的 可 行 
点 ， 提 出 了 利用 
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5 
上 
写 
十 
村 太 


利用 循环 谱 在 频谱 区 域内 具有 稀 玻 性 的 串 


压缩 感知 的 思想 对 信号 的 循环 谱 截 面 进行 信号 重 构 研 究 。 本 文 


的 分 段 正 交 匹配 追踪 算法 (STOMP) 是 在 正 交 匹配 追踪 算 


法 的 思想 基础 上 加 入 了 分 段 查询 的 思想 ， 即 在 每 次 迭代 时 设置 
一 个 阀 值 ， 来 筷 选 原子 ， 保 证 了 其 能 够 更 加 精准 地 重 构 原始 信 


1 


三) 
本 


访 真 结果 表明 在 选择 阀 值 为 2.6 情况 下 ， 重 构 信和 号 需要 的 


观测 值 较 少 ， 可 以 降低 算法 计算 量 。 在 信 噪 比 -3dB 的 情况 下 ， 
童 号 的 载 频 估计 与 实验 提供 载 频 一 致 ， 对 载波 频率 的 精度 可 以 
达到 103Hz 以 下 。 该 方法 在 工程 上 具有 一 定 的 实用 价值 。 
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